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i 1 1 Abstract 

t-H ■ 

■ Indirect information about the possible scale of supersymmetry (SUSY) breaking is provided by 

<n ; 

iy-j \ S-physics observables (BPO) as well as electroweak precision observables (EWPO). We combine the 

constraints imposed by recent measurements of the BPO BR(6 — ► 57), BR(£> S — > BR(£? U — » 

\Q ■ tv t ) and AMb s with those obtained from the experimental measurements of the EWPO Myy, 

^ ■ sin 2 6> e ff , Tz, (g — 2) M and M^, incorporating the latest theoretical calculations of these observables 

within the Standard Model and supersymmetric extensions. We perform a x 2 fit to the parameters 
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of the constrained minimal supersymmetric extension of the Standard Model (CMSSM), in which 
the SUSY-breaking parameters are universal at the GUT scale, and the non-universal Higgs model 
(NUHM), in which this constraint is relaxed for the soft SUSY-breaking contributions to the Higgs 
masses. Assuming that the lightest supersymmetric particle (LSP) provides the cold dark matter 
density preferred by WMAP and other cosmological data, we scan over the remaining parameter 
space. Within the CMSSM, we confirm the preference found previously for a relatively low SUSY- 
breaking scale, though there is some slight tension between the EWPO and the BPO. In studies 
of some specific NUHM scenarios compatible with the cold dark matter constraint we investigate 
{M Aitan. (3) planes and find preferred regions that have values of x 2 somewhat lower than in the 
CMSSM. 
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1 Introduction 



The dimensionality of the parameter space of the minimal supersymmetric extension of 
the Standard Model (MSSM) [1, 2] is so high that phenomenological analyses often make 
simplifying assumptions that reduce drastically the number of parameters. One assumption 
that is frequently employed is that (at least some of) the soft SUSY-breaking parameters 
are universal at some high input scale, before renormalization. One model based on this 
simplification is the constrained MSSM (CMSSM), in which all the soft SUSY-breaking 
scalar masses m are assumed to be universal at the GUT scale, as are the soft SUSY- 
breaking gaugino masses m-1/2 and trilinear couplings A . The assumption that squarks and 
sleptons with the same gauge quantum numbers have the same masses is motivated by the 
absence of identified supersymmetric contributions to flavour-changing neutral interactions 
and rare decays (see Ref. [3] and references therein). Universality between squarks and 
sleptons with different gauge interactions may be motivated by some GUT scenarios [4]. 
However, the universality of the soft SUSY-breaking contributions to the Higgs scalar masses 
is less motivated, and is relaxed in the non-universal Higgs model (NUHM) [5-7]. 

There are different possible approaches to analyzing the reduced parameter spaces of the 
CMSSM and the NUHM. One minimal approach would be to approximate the various theo- 
retical, phenomenological, experimental, astrophysical and cosmological constraints naively 
by 9 functions, determine the domains of the SUSY parameters allowed by their combina- 
tion, and not attempt to estimate which values of the parameters might be more or less 
likely. This approach would perhaps be adequate if one were agnostic about the existence of 
low-energy SUSY. On the other hand, if one were more positive about its existence, and keen 
to find which SUSY parameter values were more 'probable', one would make a likelihood 
analysis and take seriously any possible hints that the Standard Model (SM) might not fit 
perfectly the available data. This is the approach taken in this paper. We perform a com- 
bined x 2 analysis of electroweak precision observables (EWPO), going beyond previous such 
analyses [8,9] (see also Ref. [10]), and of 5-physics observables (BPO), including some that 
have not been included before in comprehensive analyses of the SUSY parameter space (see, 
however, Ref. [11]). In the past, the set of EWPO included in such analyses have been the 
W boson mass M^y, the effective leptonic weak mixing angle sin 2 # e ff> the anomalous mag- 
netic moment of the muon (g — 2) M , and the mass of the lightest MSSM Higgs boson mass 
Mh- Since our previous study, the theoretical link between experimental observables and 
sin 2 # c flf within the Standard Model has become more precise, changing the x 2 distribution 
for the possible MSSM contribution. We also include in this analysis a new EWPO, namely 
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the total Z boson width Y z- In addition, we now include four BPO: the branching ratios 
BR(6 — > sj), BR(B S — > fi + fi~) and BR(B U — > tu t ), and the B s mass mixing parameter 
AM Ba - For each observable, we construct the x 2 function including both theoretical and 
experimental systematic uncertainties, as well as statistical errors. The largest theoretical 
systematic uncertainty is that in BR(6 — > sj), mainly associated with the renormalization- 
scale ambiguity. Since this is not a Gaussian error, we do not add it in quadrature with the 
other errors. Instead, in order to be conservative, we prefer to add it linearly. 

For our CMSSM analysis, the fact that the cold dark matter density is known from as- 
trophysics and cosmology with an uncertainty smaller than 10 % fixes with proportional pre- 
cision one combination of the SUSY parameters, enabling us to analyze the overall x 2 value 
as a function of TOi/ 2 for fixed values of tan/? and A . The value of is fixed by the elec- 
troweak vacuum conditions, the value of m Q is fixed with a small error by the dark matter 
density, and the Higgs mass parameters are fixed by the universality assumption. As in 
previous analyses, we consider various representative values of A oc m\ji for the specific 
choices tan/3 = 10,50. Also as previously, we find a marked preference for relatively small 
values of mi/ 2 ~ 300, 600 GeV for tan/3 = 10, 50, respectively, driven largely by (g — 2) M with 
some assistance from Mw- This preference would have been more marked if the BPO were 
not taken into account. Indeed, there is a slight tension between the EWPO and the BPO, 
with the latter disfavouring smaller m^, particularly for large tan/3. As corollaries of this 
analysis, we present the x 2 distributions for the masses of various MSSM particles, includ- 
ing the lightest Higgs boson mass M h . This shows a strong preference for M h ~ 115 GeV, 
allowing M h as high as 120 GeV with A% 2 ~ 4. 

In view of the slight tension between the EWPO and BPO within the CMSSM, we have 
gone on to explore the NUHM, which effectively has Ma and /i as additional free parame- 
ters as compared to the CMSSM. In particular, we have investigated whether the NUHM 
reconciles more easily the EWPO and BPO, and specifically whether there exist NUHM 
points with significantly lower x 2 ■ As pointed out previously, generic NUHM parameter 
planes in which the other variables are held fixed do not satisfy the cold dark matter density 
constraint imposed by WMAP et al. In this paper, we introduce 'WMAP surfaces', which 
are (M^,tan/3) planes across in which the other variables are adjusted continuously so as to 
maintain the LSP density within the WMAP range. We then examine the x 2 values of the 
EWPO and BPO in the NUHM as functions over these WMAP surfaces 1 . In each of the 

1 A more complete characterization of these WMAP surfaces will be given elsewhere [12], as well as 
a discussion of their possible use as 'benchmark scenarios' for evaluating the prospects for MSSM Higgs 
phenomenology at the Tevatron, the LHC and elsewhere. 
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WMAP surfaces we find localized regions preferred by the EWPO and BPO and, in some 
cases, the minimum value of x 2 is significantly lower than along the WMAP strips in the 
CMSSM, indicating that the NUHM may help resolve the slight tension between the EWPO 
and the BPO. We explore this possibility further by investigating lines that explore further 
the NUHM parameter space in neighbourhoods of the low-% 2 points in the WMAP surfaces. 

In Sect. 2 we review the current status of the EWPO and BPO that we use, our treatment 
of the available theoretical calculations and their errors, as well as their present experimen- 
tal values. The analysis within the CMSSM can be found in Sect. 3, while the NUHM 
investigation is presented in Sect. 4. Sect. 5 summarizes our principal conclusions. 

2 Current Experimental Data 

The relevant data set includes five EWPO: the mass of the W boson, M w , the effective 
leptonic weak mixing angle, sin 2 # eff , the total Z boson width, T z , the anomalous magnetic 
moment of the muon, (g — 2) M , and the mass of the lightest MSSM Higgs boson, M h . In 
addition, we include four BPO: the branching ratios BR(6 — > S7), BR(_B S — > and 
BR(B U — > tv t ) as well as the B s mass-mixing parameter AM Bs . A detailed description of 
the EWPO and BR(6 -> s 7 ) can be found in Refs. [8,9, 13, 14]. 

In this Section we start our analysis by recalling the current precisions of the experimental 
results and the theoretical predictions for all these observables. We also display the CMSSM 
predictions for the EWPO (where new results are available), and also for the BPO. These 
predictions serve as examples of the expected ranges of the EWPO and BPO values once 
SUSY corrections are taken into account. 

In the following, we refer to the theoretical uncertainties from unknown higher-order 
corrections as 'intrinsic' theoretical uncertainties and to the uncertainties induced by the 
experimental errors of the SM input parameters as 'parametric' theoretical uncertainties. 
We do not discuss here the theoretical uncertainties in the renormalization-group running 
between the high-scale input parameters and the weak scale, see Ref. [15] for a recent dis- 
cussion in the context of calculations of the cold dark matter (CDM) density. At present, 
these uncertainties are less important than the experimental and theoretical uncertainties in 
the precision observables. 

Assuming that the nine observables listed above are uncorrelated, a x 2 fit has been 
performed with 
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71=1 




+ Xm,+Xb s - (1) 



Here i?^ xp denotes the experimental central value of the nth observable (M W) sin 2 6> e fr, T z , 
(g - 2) M and BR(6 -> s 7 ), BR(£ U -> tis t ), AM Bs ), R^ eo is the corresponding MSSM pre- 
diction and a n denotes the combined error, as specified below. Additionally, <r™ m is the 
minimum combined error over the parameter space of each data set as explained below, and 
XM h and Xb 3 denote the x 2 contribution coming from the experimental limits on the lightest 
MSSM Higgs boson mass and on BR(B S — > u + u~), respectively, which are also described 
below. 

We also list below the parametric uncertainties in the predictions on the observables 
induced by the experimental uncertainty in the top- and bottom-quark masses. These er- 
rors neglect, however, the effects of varying m t and m b on the SUSY spectrum that are 
induced via the RGE running. In order to take the m t and m b parametric uncertainties 
correctly into account, we evaluate the SUSY spectrum and the observables for each data 
point first for the nominal values m t = 171.4 GeV [16] 2 and m b (m b ) = 4.25 GeV, then for 
m t = (171.4 + 1.0) GeV and m b (m b ) = 4.25 GeV, and finally for m t = 171.4 GeV and 
m b (m b ) = (4.25 + 0.1) GeV. The latter two evaluations are used by appropriate rescal- 
ing to estimate the full parametric uncertainties induced by the experimental uncertainties 
Srrf^ = 2.1 GeV [16] 3 and 5m b (m b ) exp = 0.11 GeV. These parametric uncertainties are 
then added to the other errors (intrinsic, parametric, and experimental) of the observables 
as described in the text below. 

We preface our discussion by describing our treatment of the cosmological cold dark 
matter density, which guides our subsequent analysis of the EWPO and BPO within the 
CMSSM and NUHM. 

2.1 Cold Dark Matter Density 

Throughout this analysis, we focus our attention on parameter points that yield the cor- 
rect value of the cold dark matter density inferred from WMAP and other data, namely 
0.094 < f2cDM^ 2 < 0.129 [18]. The fact that the density is relatively well known restricts 
the SUSY parameter space to a thin, fuzzy 'WMAP hyper surf ace', effectively reducing its 



2 Using the most recent experimental value, m t = 170.9 GeV [17] would have a minor impact on our 
analysis. 

3 Using the most recent experimental m t error of 1.8 GeV [17] would also have a minor impact on our 
analysis. 
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dimensionality by one. The variations in the EWPO and BPO across this hypersurface may 
in general be neglected, so that we may treat the cold dark matter constraint effectively as 
a 5 function. For example, in the CMSSM we focus our attention on 'WMAP lines' in the 
(m 1 /2,m ) planes for discrete values of the other SUSY parameters tan/3 and A [19,20]. 
Correspondingly, in the following, for each value of mi/ 2 , we present theoretical values for 
the EWPO and BPO corresponding to the values of m on WMAP strips. 

We note, however, that for any given value of TOi/ 2 there may be more than one value 
of m that yields a cold dark matter density within the allowed range, implying that there 
may be more than one WMAP line traversing the the (m^jMo) plane. Specifically, in the 
CMSSM there is, in general, one WMAP line in the coannihilation/rapid-annihilation funnel 
region and another in the focus-point region, at higher m . Consequently, each EWPO and 
BPO may have more than one value for any given value of mi/ 2 . In the following, we restrict 
our study of the upper WMAP line to the part with m < 2000 GeV for tan/3 = 10 and 
m < 3000 GeV for tan (3 = 50, restricting in turn the range of mi/ 2 . 

The NUHM, with M A and fx, has two more parameters than the CMSSM, which charac- 
terize the degrees of non-universality of the two Higgs masses. The WMAP lines therefore 
should, in principle, be generalized to three-volumes in the higher-dimensional NUHM pa- 
rameter space where the cold dark matter density remains within the WMAP range. We 
prefer here to focus our attention on 'WMAP surfaces' that are slices through these three- 
volumes with specific fixed values for (combinations of) the other NUHM parameters. These 
WMAP surfaces are introduced in more detail in the subsequent section describing our 
NUHM analysis, and will be discussed in more detail in Ref. [12]. 

In regions that depend sensitively on the input values of m t and m^mb), such as the 
focus-point region [21] in the CMSSM, the corresponding parametric uncertainty can become 
very large. In essence, the 'WMAP hypersurface' moves significantly as m t varies (and to a 
lesser extent also m b (mb)), but remains thin. Incorporating this large parametric uncertainty 
naively in eq. (1) would artificially suppress the overall x 2 value for such points. This artificial 
suppression is avoided by adding the second term in eq. (1), where cr™ in is the value of the 
combined error evaluated for parameter choices which minimize x\ over the full data set. 

2.2 The W Boson Mass 

The W boson mass can be evaluated from 
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where a is the fine structure constant and Gf the Fermi constant. The radiative corrections 
are summarized in the quantity Ar [22]. The prediction for M w within the SM or the MSSM 
is obtained by evaluating Ar in these models and solving eq. (2) for M w . 

We use the most precise available result for M w in the MSSM [23]. Besides the full SM 
result, for the MSSM it includes the full set of one-loop contributions [23-25] as well as the 
corrections of 0(aa s ) [26] and of 0(a^ b ) [27,28] to the quantity Ap; see Ref. [23] for details. 

The remaining intrinsic theoretical uncertainty in the prediction for Mw within the 
MSSM is still significantly larger than in the SM. For realistic parameters it has been esti- 
mated as [28] 

AM intr,currcnt < 1Q ^ ^ 

depending on the mass scale of the supersymmetric particles. The parametric uncertainties 
are dominated by the experimental error of the top-quark mass and the hadronic contribu- 
tion to the shift in the fine structure constant. Their current errors induce the following 
parametric uncertainties [14] 

^current = ^ Ge y ^ AM para,m t ,currcnt _ ^ ? ^ 

cHAcC d rcnt ) = 35 x 1(T 5 AMPT' Aahad ' current ~ 6.3 MeV . (5) 
The present experimental value of Mw is [29-33] 

M exp,current = gQ ggg ± Q Q25 GeV . (6) 

We add the experimental and theoretical errors for Mw in quadrature in our analysis. 

The current status of the MSSM prediction and the experimental resolution is shown in 
Fig. 1. We note that the CMSSM predictions for Mw in the coannihilation and focus-point 
regions are quite similar, and depend little on A . We also see that small values of mi/ 2 are 
slightly preferred, reflecting the familiar fact that the experimental value of Mw is currently 
somewhat higher than the SM prediction. 

2.3 The Effective Leptonic Weak Mixing Angle 

The effective leptonic weak mixing angle at the Z boson peak can be written as 

sin 2 9 cS = \ ( 1 - Re , (7) 

where v e s and a c s denote the effective vector and axial couplings of the Z boson to charged 
leptons. We use the most precise available result for sin 2 9 e s in the MSSM [14]. The prediction 
contains the same classes of higher-order corrections as described in Sect. 2.2. 
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Figure 1: The CMSSM predictions for M\y are shown as functions of along the WMAP 
strips for (a) tan (3 = 10 and (b) tan (3 = 50 for various A values. In each panel, the centre 
(solid) line is the present central experimental value, and the (solid) outer lines show the 
current ±l-cr range. The dashed lines correspond to the full error including also parametric 
and intrinsic uncertainties. 



In the MSSM with real parameters, the remaining intrinsic theoretical uncertainty in the 
prediction for sin 2 9 e g has been estimated as [28] 

A. 2 /lintr, current ^ » i n— 5 /o\ 
sin 9 cS ' < 7 x 10 °, (8) 

depending on the SUSY mass scale. The current experimental errors of m t and Aa^ induce 
the following parametric uncertainties [14] 

^current = ^ Qe y ^ A ^2 QPpw, current _ g g x ^ (g) 

5(Aa^ cnt ) = 35 x 10~ 5 A sin 2 e p e P' Aa ^ d ' current « 12 x 10~ 5 . (10) 

The experimental value is [29,30] 

sin 2 ^p-™ 1 * = 0.23153 ± 0.00016 . (11) 

We add the experimental and theoretical errors for sin 2 6 e g in quadrature in our analysis. 

As compared with our older analyses [8,9] we now use a new result for sin 2 eff , obtained 
recently, that differs non-negligibly from that used previously, due to the inclusion of more 
higher-order corrections (which also result in a smaller intrinsic error). The corresponding 
new results in the CMSSM are shown in Fig. 2 for tan (3 = 10 (left) and tan (3 = 50 (right) as 
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functions of m\n- Whereas previously the agreement with the experimental result was best 
for mi/2 ~ 300 GeV, we now find best agreement for large m^j values. However, taking 
all uncertainties into account, the deviation for mi/2 generally stays below the level of one 
sigma. We note that the predictions for sin 2 9 e f[ in the coannihilation and focus-point regions 
are somewhat different. 
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Figure 2: The CMSSM predictions for sin 2 # e g- as functions of m 1( / 2 along the WMAP strips 
for (a) tan/3 = 10 and (b) tan/3 = 50 for various A values. In each panel, the centre (solid) 
line is the present central experimental value, and the (solid) outer lines show the current 
±l-cr range. The dashed lines correspond to the full error including also parametric and 
intrinsic uncertainties. 



2.4 The Total Z Boson Decay Width 

The total Z boson decay width, Fz, is given by 

T z = r l + r h + r^o , (12) 

where are the rates for decays into SM leptons and quarks, respectively, and T^o denotes 
the decay width to the lightest neutralino. We have checked that, for the parameters analyzed 
in this paper, always T^o = 0. However, SUSY particles enter via virtual corrections to Ti 
and r\. We use the most precise available result for Y z in the MSSM [14]. The prediction 
contains the same classes of MSSM higher-order corrections as described in Sect. 2.2. 
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So far no estimate has been made of the intrinsic uncertainty in the prediction for Yz in 
the MSSM. Following the numerical analysis in Ref. [14], we use a conservative value of 

Ar intr,current < ^ Q Me y ( 13 ) 

The current experimental errors of m t and Aa had induce the following parametric uncertain- 
ties [14] 

5m™ t = 2.1GeV A rP ara ' mt ' currcnt « 0.51 MeV, (14) 
5(Aa c d cnt ) = 35 x 10~ 5 Ar| ara ' Aahad ' currcnt « 0.32 MeV. (15) 
The experimental value is [29,30] 

r cxp,currc„t = ± 3 3 _ ^ 

We add the experimental and theoretical errors for T z in quadrature in our analysis. 

A comparison of the MSSM prediction with the experimental value is shown in Fig. 3. 
We see that the experimental value is within ~ 1/2 a standard deviation of the CMSSM 
value at large m^, which corresponds to the SM value with the same Higgs boson mass. 
The marginal improvement in the CMSSM prediction at small mi/ 2 is not significant. We 
note that the predictions for T z in the coannihilation and focus-point regions are somewhat 
different. 

2.5 The Anomalous Magnetic Moment of the Muon 

The SM prediction for the anomalous magnetic moment of the muon (see Refs. [34-38] for 
reviews) depends on the evaluation of QED contributions (see Refs. [39,40] for recent up- 
dates), the hadronic vacuum polarization and light-by-light (LBL) contributions. The former 
have been evaluated in Refs. [38,41-46] and the latter in Refs. [47-51]. The evaluations of 
the hadronic vacuum polarization contributions using e + e~ and r decay data give somewhat 
different results. In view of the fact that recent e + e~ measurements tend to confirm ear- 
lier results, whereas the correspondence between previous r data and preliminary data from 
BELLE is not so clear, and also in view of the additional uncertainties associated with the 
isospin transformation from r decay, we use here the latest estimate based on e + e~ data [46]: 

a^ eo = (11 659 180.5 ± 4.4 had ± 3.5 LBL ± 0.2 QED+EW ) x lO" 10 , (17) 

where the source of each error is labeled. We note that the new e + e~ data sets that have 
recently been published in Refs. [52-54] have been partially included in the updated estimate 

of (g - 2),. 
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Figure 3: The CMSSM predictions for V z as functions of m^j along the WMAP strips for 
(a) tan j3 = 10 and (b) tan /3 = 50 for various A values. In each panel, the centre (solid) line 
is the present central experimental value, and the (solid) outer lines show the current ±l-cr 
range. The dashed lines correspond to the full error including also parametric and intrinsic 
uncertainties. 



The SM prediction is to be compared with the final result of the Brookhaven (g — 2) M 
experiment E821 [55,56], namely: 

a™? = (11 659 208.0 ± 6.3) x HT 10 , (18) 

leading to an estimated discrepancy [46,57] 

G cx P _ a theo = ( 27>5 ± g ^ x 10 -io ? ( 19 ) 

equivalent to a 3.3-u effect 4 . While it would be premature to regard this deviation as a firm 
evidence for new physics, within the context of SUSY, it does indicate a preference for a 
non-zero contribution. 

Concerning the MSSM contribution, the complete one-loop result was evaluated a decade 
ago [58]. In view of the correlation between the signs of (g — 2) M and of [i [59], variants of the 
MSSM with fi < are already severely challenged by the present data on a M , whether one 
uses either the e + e~ or r decay data, so we restrict our attention in this paper to models with 
fi > 0. In addition to the full one-loop contributions, the leading QED two-loop corrections 

4 Three other recent evaluations yield slightly different numbers [37,38,43], but similar discrepancies with 
the SM prediction. 



10 



have also been evaluated [60]. Further corrections at the two-loop level have been obtained 
recently [61,62], leading to corrections to the one-loop result that are < 10%. These 
corrections are taken into account in our analysis according to the approximate formulae 
given in Refs. [61,62]. 

The current status of the CMSSM prediction and the experimental resolution is shown 
in Fig. 4, where the 1- and 2-cr bands are shown. We note that the coannihilation and focus- 
point region predictions for a M are quite different. For tan/? = 10, the focus-point prediction 
agrees less well with the data, whereas for tan (3 = 50 the focus-point prediction does agree 
well in a limited range of ~ 200 GeV. 
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Figure 4: The CMSSM predictions for (g — 2) M , Aa M , as functions of along the WMAP 
strips for (a) tan (3 — 10 and (b) tan (3 = 50 for various A values. In each panel, the centre 
(solid) line is the present central experimental value, and the solid (dotted) outer lines show 
the current ±1(2)-<t ranges. 



2.6 The Mass of the Lightest MSSM Higgs Boson 

The mass of the lightest CP-even MSSM Higgs boson can be predicted in terms of the other 
MSSM parameters. At the tree level, the two CP-even Higgs boson masses are obtained as 
functions of Mz, the CP-odd Higgs boson mass M^, and tan/?, whereas other parameters 
enter into the loop corrections. We employ the Feynman-diagrammatic method for the the- 
oretical prediction of Mh, using the code FeynHiggs [63-65], which includes all numerically 
relevant known higher-order corrections. The status of these results can be summarized 
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as follows. For the one-loop part, the complete result within the MSSM is known [66-68]. 
Computation of the two-loop effects is quite advanced: see Ref. [69] and references therein. 
These include the strong corrections at 0(a t a s ) and Yukawa corrections at 0{af) to the 
dominant one-loop 0(a t ) term, and the strong corrections from the bottom/sbottom sec- 
tor at 0{ai ) a s ). In the case of the b/b sector corrections, an all-order resummation of the 
tan (3 -enhanced terms, 0(otb(a s tan/5) n ), is also known [70,71]. Most recently, the 0(a t ab) 
and 0(al) corrections have been derived [72] 5 . The current intrinsic error of M h due to 
unknown higher-order corrections has been estimated to be [13,69,75,76] 

AM intr,current = 3 ? (2Q) 

which we interpret effectively as a ~ 95 % confidence level limit: see below. 

It should be noted that, for the unconstrained MSSM with small values of M A and values 
of tan/3 which are not too small, a significant suppression of the hZZ coupling can occur 
compared to the SM value, in which case the experimental lower bound on Mh may be more 
than 20 GeV below the SM value [77]. However, we have checked that within the CMSSM 
and the other models studied in this paper, the hZZ coupling is always very close to the 
SM value. Accordingly, the bounds from the SM Higgs search at LEP [78] can be taken over 
directly (see e.g. Refs. [79,80]). 

Concerning the x 2 analysis, we use the complete likelihood information available from 
LEP. Accordingly, we evaluate as follows the M h contribution to the overall x 2 function 6 . 
Our starting points are the CL s (M h ) values provided by the final LEP results on the SM 
Higgs boson search, see Fig. 9 in Ref. [78] 7 . We obtain by inversion from CL s (M h ) the 
corresponding value of x 2 (M h ) determined from Ref. [81] 

ierfc(^/|xW0) = CL s (M h ) , (21) 

and note the fact that CL s (M h = 116.4 GeV) = 0.5 implies that x 2 (116.4 GeV) = as 
is appropriate for a one-sided limit. Correspondingly we set x 2 (M ft > 116.4 GeV) = 0. 
The theoretical uncertainty is included by convolving the likelihood function associated with 
X 2 (M/i) and a Gaussian function, $i. 5 (x), normalized to unity and centered around Mh, 
whose width is 1.5 GeV: 

X 2 (M h ) = -2 log e-* 2 ^ 2 $ 1 . 5 (M, - x) dx) . (22) 

5 A two- loop effective potential calculation has been presented in Ref. [73], including now even the leading 
threedoop corrections [74], but no public code based on this result is currently available. 
6 We thank P. Bechtle and K. Desch for detailed discussions and explanations. 
7 We thank A. Read for providing us with the CL S values. 
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In this way, a theoretical uncertainty of up to 3 GeV is assigned for ~ 95% of all Mh values 
corresponding to one parameter point. The final x\i h is then obtained as 

XM h = x\M h )-x\llQAGeY) for M h < HQ A GeV , (23) 

X 2 Mh = for M h > 116.4 GeV , (24) 

and is then combined with the corresponding quantities for the other observables we consider, 
see eq. (1). 

We show in Fig. 5 the predictions for Mh in the CMSSM for tan/3 = 10 (left) and 
tan/3 = 50 (right). The predicted values of Mh are similar in the coannihilation and focus- 
point regions. They depend significantly on A , particularly in the coannihilation region, 
where negative values of A tend to predict very low values of that are disfavoured by 
the LEP direct search. Also shown in Fig. 5 is the present nominal 95 % C.L. exclusion limit 
for a SM-like Higgs boson, namely 114.4 GeV [78], and a hypothetical LHC measurement 
of Mh = 116.4 ± 0.2 GeV. We recall that we use the numerical value of the LEP Higgs 
likelihood function in our combined analysis. 
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Figure 5: The CMSSM predictions for Mh as functions of mi/2 with (a) tan/3 = 10 and (b) 
tan/? = 50 for various Aq. We also show the present 95% C.L. exclusion limit of 114.4 GeV 
and a hypothetical LHC measurement of Mh = 116.4 ± 0.2 GeV. 



2.7 The decay b — > S7 

Since this decay occurs at the loop level in the SM, the MSSM contribution might a priori be 
of similar magnitude. A recent theoretical estimate of the SM contribution to the branching 
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ratio at the NNLO QCD level is [82] 



BR(b -> s 7 ) = (3.15 ± 0.23) x 10 



-4 



(25) 



We record that the error estimate for BR(6 — > S7) is still under debate, and that other SM 
contributions to b — > S7 have been calculated Refs. [83,84], but these corrections are small 
compared with the theoretical uncertainty quoted in (25). 



For comparison, the present experimental value estimated by the Heavy Flavour Aver- 
aging Group (HFAG) is [3,85] 



where the first error is the combined statistical and uncorrelated systematic uncertainty, the 
latter two errors are correlated systematic theoretical uncertainties and corrections respec- 
tively. 

Our numerical results have been derived with the BR(6 — > S7) evaluation provided in 
Refs. [86-88], incorporating also the latest SM corrections provided in Ref. [82]. The calcula- 
tion has been checked against other approaches [89-91]. For the current theoretical intrinsic 
uncertainty of the MSSM prediction for BR(6 — > 57) we use the SM uncertainty given in 
eq. (25) and add linearly the intrinsic MSSM corrections 0.15 x 10~ 4 [89,91] and the last two 
errors given by HFAG of ~ 0.13 x 10~ 4 [3]. The full intrinsic error is then added linearly to 
the sum in quadrature of the experimental error given by HFAG as 0.24 and the parametric 
error. 

In Fig. 6 we show the predictions in the CMSSM for BR(6 — > 57) for tan/? = 10,50 
as functions of m^, compared with the 1-a experimental error (full line) and the full error 
(dashed line, but assuming a negligible parametric error). For tan (3 = 10, we see that positive 
values of A are disfavoured at small my 2, an d that small values of my 2 are disfavoured for 
all the studied values of A if tan (3 = 50. 

2.8 The Branching Ratio for B s — > 

The SM prediction for this branching ratio is (3.4 ± 0.5) x 10 -9 [92], and the present exper- 
imental upper limit from the Fermilab Tevatron collider is 1.0 x 10~ 7 at the 95% C.L. [93], 
providing ample room for the MSSM to dominate the SM contribution. The current Teva- 
tron sensitivity is based on an integrated luminosity of about 780 pb -1 collected at CDF. The 



BR(b -> s 7 ) = (3.55 ± 0.24t° f ± 0.03) x 10~ 4 , 



(26) 
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Figure 6: The CMSSM predictions for BR(6 — > 57) as functions of raj/2 along the WMAP 
strips for tan/3 = 10 (left) and tan/3 = 50 (right) and various choices of Aq. The central 
(solid) line indicates the current experimental central value, and the other solid lines show 
the current ±1-<t experimental range. The dashed line is the ±l-cr error including also the 
full intrinsic error (see text). 



exclusion bounds can be translated into a x 2 function for each value of BK(B S 

X 2 (B S ) = x 2 (BR(£? s - ^V)) , 



8. 



(27) 



with x 2 (BR(i? s — > fi + jj,—) < 0.266 x 10 -7 ) = 0. The theory uncertainty is included by 
convolving the likelihood function associated with x 2 {B s ) and a Gaussian function, $ t h(^), 
normalized to unity and centered around BR(13 S — > /i + /i _ ), whose width is given by the 
theory uncertainty, see below. Consequently, 



X \B S ) = -2 log ( / e-* 2 M/ 2 $ th (BR(B s - - *) dx 



The final Xb s is then obtained as 

Xl 3 = X 2 (B S ) ~ X 2 (0.266 x 10- 7 ) for BR(B S 

X 2 Bs =0 for BR(B S 



> 0.266 x 10~ 7 
12,+fT) < 0.266 x 10~ 7 



(28) 



(29) 
(30) 



The Tevatron sensitivity is expected to improve significantly in the future. The limit that 
could be reached at the end of Run II is ~ 2 x 10~ 8 assuming 8 fb _1 collected with each 



8 We thank C.-J. Stephen and M. Herndon for providing the x 2 numbers. A slightly more stringent 
upper limit of 0.93 x 10~ 7 at the 95% C.L. has been announced more recently by the DO Collaboration [94]. 
However, the corresponding x 2 function is not available to us. Since the difference to the result employed 
here is small, we expect only a minor impact on our analysis. 
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detector [95]. A sensitivity even down to the SM value can be expected at the LHC. Assuming 
the SM value, i.e. BR(5 S — > ss 3.4 x 10~ 9 , it has been estimated [96] that LHCb can 

observe 33 signal events over 10 background events within 3 years of low-luminosity running. 
Therefore this process offers good prospects for probing the MSSM. 

For the theoretical prediction we use results from Ref. [97], which are in good agreement 
with Ref. [98] . This calculation includes the full one- loop evaluation and the leading two-loop 
QCD corrections. 

The theory error is estimated as follows. We take into account the parametric uncertainty 
induced by [99] 

f Bs = 230 ± 30 MeV . (31) 
The most important SUSY contribution to BR(£> S — > scales as 

BR(B S -> /iV) ~ J| • (32) 

In the models that predict the value of Ma at the low-energy scale, i.e. in our case the 
CMSSM, we additionally include the parametric uncertainty due to the shift in Ma in eq. (32) 
that is induced by the experimental errors of m t and in the RGE running [98]. These 
errors are added in quadrature. The intrinsic error is estimated to be negligible as compared 
to the parametric error. Thus the parametric error constitutes our theory error entering in 
eq. (28). 

In Fig. 7 the CMSSM predictions for BK(B S — > /U + /i~) for tan (3 = 10,50 as functions of 
TO1/2 are compared with the present Tevatron limit. For tan/5 = 10 (left plot) the CMSSM 
prediction is significantly below the present and future Tevatron sensitivity. However, already 
with the current sensitivity, the Tevatron starts to probe the CMSSM coannihilation region 
for tan/9 = 50 and A > 0, whereas the CMSSM prediction in the focus-point region is 
significantly below the current sensitivity. 

2.9 The Branching Ratio for B u — > rv T 

The decay B u — > rv T has recently been observed by BELLE [100], and the experimental 
world average is given by [11, 100, 101] 

BR(B U -> rz/ r ) cxp = (1.31 ± 0.49) x 10~ 4 . (33) 

We follow Ref. [102] for the theoretical evaluation of this decay. The main new contribution 
within the MSSM comes from the direct-exchange of a virtual charged Higgs boson decaying 
into ru T . Taking into account the resummation of the leading tan/5 enhanced corrections, 
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Figure 7: The CMSSM predictions for BR(B S — ► as functions of along the 

WMAP strips for tan/3 = 10 (left) and tan/3 = 50 (right) and various choices of Aq. The 
solid line indicates the current experimental 95% CL. exclusion bound. 



within scenarios with minimal flavor violation such as the CMSSM and the NUHM, the ratio 
of the MSSM result over the SM result can be written as 



BR(B U — » tz/ t )mssm 
BR(B U -> tu t ) sm 



m 2 B \ tan 2 /3 



(34) 



\M 2 H± ) 1 + e tan/3_ 

Here £o denotes the effective coupling of the charged Higgs boson to up- and down-type 
quarks, see Ref. [102] for details. The deviation of the experimental result from the SM 
prediction can be expressed as 

||^^% = 0.93 ±0.41, (35) 
BR{B U -> Tv r ) SM 

where the error includes the experimental error as well as the parametric errors from the 
various SM inputs. We use eq. (34) for our theory evaluation, which can then be compared 
with eq. (35), provided that the value of AM# d agrees sufficiently well in the SM and in the 
MSSM (which we assume here). As an error estimate we use the combined experimental and 
parametric error from eq. (35), an estimated intrinsic error of ~ 2%, and in the CMSSM, 
as for BR(S S — > /i + /i - ), an additional parametric error from Mh±, evaluated from RGE 
running. These errors have been added in quadrature. 

We show in Fig. 8 the theoretical results for the ratio of CMSSM/SM for BR(5 U — ► tu t ) 
as functions of for tan/3 = 10,50. These results are also compared with the present 
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experimental result. The central (solid) line indicates the current experimental central value, 
and the other solid (dotted) lines show the current ±l(2)-cx ranges from eq. (35). For tan f3 = 
10 the SM result is reproduced over most of the parameter space. Only very small mi/2 values 
give a ratio visibly smaller than 1. For tan /3 = 50 the result varies strongly between and 1, 
and the CMSSM could easily account for the small deviation of the central value of the 
experimental result from the SM prediction, should that become necessary. The prediction 
in the focus-point region is somewhat closer to the SM value. 
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Figure 8: The predictions for the ratio CMSSM/SM for BK(B U — > ru T ) as functions of m 1( / 2 
along the WMAP strips for tan (3 — 10 (left) and tan (3 = 50 (right) and various choices of 
A . The central (solid) line indicates the current experimental central value, and the other 
solid (dotted) lines show the current ±1(2)-<t ranges. 



2.10 The B s -B s Mass Difference AM Bs 

The B s -B s oscillation frequency and consequently the the B s -B s mass difference has recently 
been measured by the CDF Collaboration [103], 

(AMeJexp = 17.77 ± 0.12 ps- 1 , (36) 

which is compatible with the broader range of the result from DO [104]. 

We follow Ref. [102] for the theory evaluation. The main MSSM contribution to the 
B s -B s oscillation comes from the exchange of neutral Higgs bosons, but we use here the full 
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result given in Ref. [102] (taken from Ref. [105]), where the leading dependence is given as 

(AMbJmssm m b (m b )m s (m h ) 



1 



(AMbJsm 

The SM value, obtained from a global fit, is given by [106] 

(AM b Jsm = 19.0 ±2.4 ps- 1 



resulting in 



(AM 



B s ) cxp 



0.93 ±0.13 . 



(37) 



(38) 



(39) 



(AMbJsm 

The error in eq. (39) is supplemented by the parametric errors in eq. (37) from m s (mb) = 93± 
17 MeV and, in the case of the CMSSM, as for BR(5 S — > an additional parametric 

error from M4. These errors are added in quadrature. The intrinsic error, in comparison, is 
assumed to be negligible. 
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Figure 9: The predictions for the ratio CMSSM/SM for AM# s as functions of m 1 / 2 along 
the WMAP strips for tan/? = 10 (left) and tan/? = 50 (right) and various choices of A . 
The central (solid) line indicates the current experimental central value, and the other solid 
(dotted) lines show the current ±1(2)-<t ranges. 

In Fig. 9 we show the results for the ratio of CMSSM/SM for AM# s as functions of 
m 1( / 2 for tan/3 = 10,50. These are also compared with the present experimental result. 
The central (solid) line indicates the current experimental central value, and the other solid 
(dotted) lines show the current ±l(2)-a ranges from eq. (39). For tan/3 = 10 the SM result 
is reproduced over the whole parameter space. Only for tan/3 = 50 and rrii/ 2 ^ 500 GeV in 
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the coannihilation region can the CMSSM prediction be significantly lower than 1. Here the 
CMSSM could account for the small deviation of the experimental result from the central 
value SM prediction, should that be necessary. 

3 CMSSM Analysis Including EWPO and BPO 

We now use the analyses of the previous Section to estimate the combined x 2 function for 
the CMSSM as a function of m^, using the master formula (1). As a first step, Fig. 10 
displays the x 2 distribution for the EWPO alone. 

In the case tan (3 = 10 (left panel of Fig. 10), we see a well-defined minimum of x 2 f° r 
mi/2 ~ 300 GeV when A > 0, which disappears for large negative A and is not present in 
the focus-point region. The rise at small mi/ 2 is due both to the lower limit on Mh coming 
from the direct search at LEP and to (g — 2) M , whilst the rise at large mi/ 2 is mainly due to 
(g — 2) M (see Fig. 4). The measurement of M w (see Fig. 1) leads to a slightly lower minimal 
value of x 2 -, but there are no substantial contributions from any of the other EWPO. The 
preference for A > in the coannihilation region is due to M h (see Fig. 5), and the relative 
disfavour for the focus-point regions is due to its mismatch with (g — 2) M (see Fig. 4). 

In the case tan/5 = 50 (right panel of Fig. 10), we again see a well-defined minimum of 
X 2 , this time for mi/ 2 ~ 400 to 500 GeV, which is similar for all the studied values of A . In 
this case, there is also a similar minimum of x 2 for the focus-point region at mi/ 2 ~ 200 GeV. 
The increase in x 2 at small rai/ 2 is due to (g — 2) M as well as M h) whereas the increase at 
large mi/ 2 is essentially due to (g — 2) M . We note that the overall minimum of % 2 ~ 2 is 
similar for both values of tan /3, and represents an excellent fit in each case. 

Fig. 11 shows the corresponding combined x 2 fc> r the BPO alone. For both values of tan (3, 
these prefer large values of m^, reflecting the fact that there is no hint of any deviation 
from the SM, and the overall quality of the fit is good. Small values of mi/ 2 are disfavoured, 
particularly in the coannihilation region with A > 0, mainly due to b — > sj. The focus- 
point region is generally in very good agreement with the BPO data, except at very low 
mi/a < 400 GeV for tan/? = 50. 

Finally, we show in Fig. 12 the combined x 2 values for the EWPO and BPO, computed 
in accordance with eq. (1). We see that the global minimum of x 2 ~ 4.5 for both values 
of tan p. This is quite a good fit for the number of experimental observables being fitted, 
and the x 2 /d.o.f. is similar to the one for the EWPO alone. This increase in the total x 2 
reflects the fact that the BPO exhibit no tendency to reinforce the preference of the EWPO 
for small m^: rather the reverse, in fact. For both values of tan/5, the focus-point region is 
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Figure 10: The combined x 2 function for the electroweak observables Mw, sin 2 ^, Tz, 
(g - 2) /t and M h , evaluated in the CMSSM for tan/5 = 10 (left) and tan/3 = 50 (right) for 
various discrete values of A . We use m t = 171.4 ±2.1 GeV and mb(mb) = 4.25 ±0.11 GeV, 
and mo is chosen to yield the central value of the cold dark matter density indicated by 
WMAP and other observations for the central values of m t and m^mb). 
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Figure 11: The combined \ 2 function for the h physics observables BR(6 — > 57), BR(5 S — > 
BR(B U -> tu t ) and AM Bs , evaluated in the CMSSM for tan/5 = 10 (left) and 
tan/? = 50 (right) for various discrete values of Aq. We use m t = 171 A ± 2.1 GeV and 
Tribimb) = 4.25 ±0.11 GeV, and mo is chosen to yield the central value of the cold dark 
matter density indicated by WMAP and other observations for the central values of m t and 
m b {m b ). 
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disfavoured by comparison with the coannihilation region, though this effect is less important 
for tan/3 = 50. For tan/3 = 10, m\j% ~ 300 GeV and A > are preferred, whereas, for 
tan /3 = 50, mi/2 ~ 600 GeV and Aq < are preferred. This change-over is largely due to the 
impact of the LEP constraint for tan /? = 10 and the b — > 37 constraint for tan j3 = 50. 
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Figure 12: The combined x 2 function for the electroweak observables M\ 
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(g — 2) M , Mh, and the b physics observables BR(6 — > 37), BR(B S — > /i + yU~), BR(_B U — > rz/ r ) 
and AMb s , evaluated in the CMSSM for tan/3 = 10 (left) and tan/3 = 50 (right) for various 
discrete values of A . We use m t = 171.4 ±2.1 GeV and mb(m b ) = 4.25 ±0.11 GeV, and 
mo is chosen to yield the central value of the cold dark matter density indicated by WMAP 
and other observations for the central values of m t and mJmi,). 



We display in Fig. 13 the \ 2 functions for various SUSY masses in the CMSSM for 
tan/3 = 10, including (a) m^o, (b) m^o and m-± (which are very similar), (c) m^, (d) Ma, 
(e) mfj and (f) nig. We see two distinct populations of points, corresponding to the x ~ ^1 
coannihilation (which is favoured) and focus-point regions (which is disfavoured). In the 
latter region, very low values of m^j are preferred, as can be seen in panels (a) and (f), 
relatively small values of /z, as can be seen in panel (b), large values of mo, as can be seen 
in panels (c) and (e), and large values of Ma, as can (not) be seen in panel (d). Compared 
to the analysis in Ref. [9], where BR(6 — > 37) was the only BPO included, and where a top 
quark mass of 172.7 GeV was used, there is no significant shift of the values of the masses 
where x 2 has its minimum, which is in the coannihilation region. As before, the present 
analysis gives hope for seeing squarks and gluinos in the early days of the LHC (panels (e) 
and (f)), and also hope for seeing charginos, neutralinos and staus at the ILC (panels (a), (b) 
and (c)), whereas observing the heavier Higgs bosons would be more challenging (panel (d)). 
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Figure 13: Various SUSY masses are presented with their respective \ 2 value in the CMSSM 
for tan/5 = 10. The panels show (a) rn^o, (b) m^o and m~± (which are very similar), (c) m^, 
(d) Ma, (e) and (f) m§. 
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In Fig. 14 we show the analogous x 2 functions for various SUSY masses in the CMSSM for 
tan/3 = 50: (a) m^o, (b) m^o and m-± (which are very similar), (c) m^, (d) Ma, (e) and 
(f) mj. We again see the clear separation between the focus-point and coannihilation regions, 
interpolated by a light-Higgs pole strip, and that the coannihilation region is somewhat 
preferred. As for lower tan/5, small values of mi/ 2 and larger values of m are preferred, 
and also small values of /j, and larger values of Ma- Again as for tan (3 = 10, compared to 
the analysis in Ref. [9], where BR(6 — > 57) was the only BPO included and where a top 
quark mass of 172.7 GeV was used, we do not find a significant shift in the values of the 
masses with lowest x 2 ■ The sparticle masses are generally higher than for tan (3 = 10: finding 
squarks and gluinos should still be 'easy' at the LHC, but seeing charginos, neutralinos and 
staus at the ILC would be more challenging, depending on its center-of-mass energy. 

Analogously to the sparticle masses in Figs. 13 and 14, we display in Fig. 15 the total 
X 2 functions for M h) as calculated in the CMSSM for tan/5 = 10 (left panel) and tan/5 = 
50 (right panel). We recall that this theoretical prediction has an intrinsic uncertainty 
of ~ ±1.5 GeV, which should be combined with the experimental error in m t . It is a clear 
prediction of this analysis that M h should be very close to the LEP lower limit, and probably 
< 120 GeV, though a value as large as ~ 123 GeV is possible (but is x 2 disfavoured), 
particularly if tan (3 = 50. 

In the case of the SM, it is well known that tension between the lower limit on M h from 
the LEP direct search and the relatively low value of M h preferred by the EWPO has recently 
been increasing [30,31]. This tension is strongly reduced within the CMSSM, particularly 
for tan/5 = 50. We display in Fig. 16 the global x 2 functions for the EWPO and BPO, but 
this time omitting the contribution for the LEP Higgs search. This corresponds to the fitted 
value of M h in the CMSSM. Comparing Fig. 16 and Fig. 15, we see that all data (excluding 
M h ) favour a value of M h ~ 110 GeV if tan/3 = 10 and M h ~ 115 GeV if tan/5 = 50. On 
the other hand, the currently best-fit value of is 76 GeV [30], i.e. substantially below 
the SM LEP bound of 114.4 GeV [78]. In comparison to the favoured values including the 
LEP limits we get a ~ 5 GeV smaller value of M h if tan/5 = 10, whereas the difference is 
only ~ 1 GeV if tan/5 = 50. 9 Correspondingly, comparing Fig. 16 and Fig. 15, we see that 
the LEP limit increases the value of x 2 by ~ 3.5 for the tan/5 = 10 case, but by only ~ 1 for 
the tan (3 = 50 case. However, we emphasize that for both cases there are quite good global 
fits to all the EWPO and BPO with X 2 ~ 4.5. 

9 We also recall that the estimated theoretical uncertainty in Mh for fixed values of the CMSSM parameters 
is ±3 GeV. 
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Figure 14: Various SUSY masses are presented with their respective \ 2 value in the CMSSM 
for tan/5 = 50. The panels show (a) rn^o, (b) m^o and m~± (which are very similar), (c) m^, 
(d) Ma, (e) and (f) m§. 



25 



10 - 



95 




CMSSM, n > 

tanp = 10, A 

tanp = 10, A 

tan|3 = 10, A 

tanp = 10, A 

tanp = 10, A 

i i i I i i 
100 



0, m, = 171.4 GeV 
= 

= +m 1/2 
= -m„ 2 
= +2 m, /2 
= -2m,. 



I 



'1/2 



105 



M, 



110 

, [GeV] 



115 



120 



14 



12 



10 - 



CO 
"O 

o 



CMSSM, |j, > 0, m, = 1 71 .4 GeV \ 
m tanp = 50, A = 

tanp = 50, A = +m 1/2 

tanp = 50, A = -m 1/2 
• tanp = 50, A = +2 m yz 
o tanp = 50, A = -2 m 1/2 




125 



95 



100 



105 



110 115 

M h [GeV] 



120 125 



Figure 15: The combined x 2 function for Mh, as obtained from the combined analysis of all 
EWPO and BPO, evaluated in the CMSSM for tan/3 = 10 (left) and tan/? = 50 (right) for 
various discrete values of A . We use m t = 171.4 ±2.1 GeV and m b {m b ) = 4.25 ±0.11 GeV, 
and mo is chosen to yield the central value of the cold dark matter density indicated by 
WMAP and other observations for the central values of mt and m^mj). 
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Figure 16: The combined x 2 function for Mh, as obtained from a combined analysis of all 
EWPO and BPO except the LEP Higgs search, as evaluated in the CMSSM for tan (5 = 10 
(left) and tan (3 = 50 (right) for various discrete values of Aq. We use m t = 171.4 ±2.1 GeV 
and mb(mb) = 4.25 ± 0.11 GeV, and mo is chosen to yield the central value of the cold dark 
matter density indicated by WMAP and other observations for the central values of m t and 
m b (m b ). 
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4 NUHM Analysis Including EWPO and BPO 



The CMSSM is a very particular case of the general MSSM. It has a manageable parameter 
space, but may not be able to capture all the possibilities available in the general MSSM. 
Specifically, as we have seen, while providing a better fit to the EWPO than the SM, it 
provides no improvement for the BPO, and there is a slighttension between the EWPO and 
the BPO within the restrictive CMSSM framework. For these and other reasons, we now 
consider the NUHM. The dimensionality of the MSSM parameter space is increased by a 
manageable amount compared to the CMSSM, namely two extra dimensions. Moreover, 
there was no strong phenomenological motivation for assuming universality for the Higgs 
masses, and there is reason to hope that relaxing the Higgs universality assumption may 
help reconcile the EWPO and BPO. As we have seen, the EWPO prefer a specific range of 
m \/2 (~ 300 GeV for tan /3 — 10 and ~ 500 GeV for tan f3 = 50) in the coannihilation region, 
and disfavour the focus-point region (particularly for tan/3 = 10). Within the CMSSM, the 
electroweak vacuum conditions fix the corresponding values of \fj,\ and Ma- These values are 
not very crucial for the EWPO, but are potentially important for the BPO. For example, 
the extra MSSM contribution to b — > 57 is small only if the supplementary charged-Higgs 
and chargino diagrams cancel to some extent, which imposes a specific condition on their 
masses. This may not be satisfied within the CMSSM, but is not incompatible a priori with 
the NUHM, in which and M4 become (to some extent) free parameters. 

Just as we focused attention in the previous CMSSM analysis on WMAP lines in pa- 
rameter space, where the cold dark matter density falls within the range allowed by WMAP 
and other astrophysical and cosmological observations, we also focus on 'WMAP surfaces' 
in in the NUHM parameter space. Many NUHM parameter planes have been considered in 
the past [5,6,107] and, as in the CMSSM, generically the dark matter constraint is satisfied 
only in thin strips in each NUHM plane. Many phenomenological studies of MSSM Higgs 
physics have analyzed the possibilities in (M^tan/?) planes under different hypotheses for 
other MSSM parameters. In particular, in the general MSSM framework, (M^,tan/3) planes 
have been suggested for phenomenological Higgs physics analyses [108-110], neglecting the 
constraints coming from CDM. In the NUHM, in order to keep the dark matter density 
within the WMAP range across generic regions of such a (Ma, tan (3) plane, one must adjust 
one or more of the free parameters continuously across the plane. We propose here two 
strategies for specifying suitable parameter scans. 

We consider first a typical (Ma, tan (3) plane for fixed mo, mi/2 and /i, as shown for 
example in Fig. 5a of Ref. [5], where the choices m = 800 GeV, mi/ 2 = 600 GeV and 
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\x = 1000 GeV were made. In this example, the relic density exceeds the WMAP upper 
limit almost everywhere in the (Ma, tan (5) plane, except along a narrow vertical strip where 
m^o ~ Ma/2, and rapid direct-channel annihilation suppresses the relic density below the 
WMAP range. On either side of this strip, there are thin regions where the relationship 
between M A and mi/ 2 is such that the relic density falls within the WMAP range whatever 
the value of tan/5. Building on this observation, we study a (M^,tan/3) plane PI with the 
same values of m = 800 GeV and /x = 1000 GeV, but with chosen to vary across the 
plane so as to maintain the WMAP relationship with M A - 

-M A - 12.5 GeV < m 1/2 < -M A + 37.5 GeV. (40) 
8 8 

As we saw earlier, within the CMSSM, smaller values of mo are preferred. We therefore con- 
sider also a (M A , tan (3) plane P2 with the fixed values m = 300 GeV and /i = 800 GeV, and 
mi/2 again adjusted continuously across the plane so as to maintain the WMAP relationship 
with Ma'- 

1.2M a - 40 GeV < m 1/2 < 1.2M A + 40 GeV. (41) 

Many more examples could be chosen, but these serve as representative examples of NUHM 
planes that enable us to explore the possibility of reducing the tension between the EWPO 
and BPO. 

We also consider two more examples, inspired by the (fx, Ma) planes also shown in Ref. [5]. 
There we see that, for fixed values of mi/ 2 and m , there is a 'magic' value of fi > which 
provides a suitable value of the relic density for almost all values of Ma, the exception being 
a narrow strip around the rapid-annihilation funnel where m^o ~ Ma/2. The value of fi 
varies with tan/3, but a suitable scan yields a (Ma, tan /3) plane where the relic density 
falls within the WMAP range for all except a sliver of M A values that broadens somewhat 
as tan/3 increases. Within this sliver, the relic density falls below the WMAP range: this 
region is therefore not incompatible with cosmology, but would require a supplementary 
source of cold dark matter. One example of such a plane, P3, has fixed mi/ 2 = 500 GeV 
and m = 1000 GeV, with \i in the range 

H = 250 - 400 GeV. (42) 

The other example P4 has fixed mi/ 2 = 300 GeV and m = 300 GeV, with /i in the range 

H = 200- 350 GeV. (43) 

The four scenarios are summarized in Tab. 1. In the analyses below, we quote the mini- 
mal values of x 2 f° r values of /i within the ranges (42), (43), for the planes P3 and P4, 
respectively. 
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Table 1: The four NUHM scenarios, with M A and tan/3 kept as free parameters. All masses 
are in GeV. 

We now consider the most important contributions to the likelihood functions for these 
four (M A) tan (3) planes. 

The principal contributions to the overall x 2 value, and hence also to the likelihood 
function, for the (M A ,tan(3) plane for scenario PI are shown in Fig. 17. We see in panel (a) 
that the \ 2 value of a M is not very satisfactory anywhere in the plane, but particularly not 
at large M A and small tan/3. Panel (b) shows that the LEP lower limit on M h disfavours 
Ma < 300 GeV in this scenario. Small values of M A are also disfavoured by b — > 57, as seen 
in panel (c), and large values of M A and tan/3 are also disfavoured, but to a lesser extent. 
Panels (d) and (e) show that large values of tan (3 and small values of M A are disfavoured 
by both B s — > and B u — > tu t . Finally, panel (f) shows the values of the combined 

EWPO and BPO x 2 function for scenario PI throughout the (M^,tan/3) plane. We see 
that the best-fit point has M A ~ 440 GeV and tan/3 ~ 50. This is the optimal compromise 
between and b — > S7 that also respects the M h and other BPO constraints. We note that 
this best-fit point has x 2 — 7.1, which is not a significant improvement, but even slightly 
worse than the CMSSM fits discussed in the previous Section. We also display in panel (f) 
of Fig. 17 the A% 2 = 2.30 and 4.61 contours, which would correspond to the 68 % and 95 % 
C.L. contours in the (M A ,teai(3) plane if the overall likelihood distribution, £ oc e - * 2 / 2 , 
was Gaussian. This is clearly only roughly the case in this analysis, but these contours 
nevertheless give interesting indications on the preferred region in the (M A ,teai(3) plane. 

We do not show results in the upper right corners of these planes (with high M A and 
high tan/3) because there the relic density in this region is low compared to the preferred 
WMAP value. In this region, for the choice of m and /i and the range of my 2 given in 
eq. (40), we are sitting too close to the funnel. However, these points could be brought into 
agreement with WMAP, by extending the sampled range in my 2 to lower values. The lower 
left portions of these planes are missing because of the finite resolution of our scan. In these 
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regions of fixed low values of M A and tan/5, the relic density is very sensitive to mi/ 2 , and 
viable points are missed with the resolution in my 2 of 10 GeV that we use. One final remark 
on Fig. 17 concerns high values of tan/3. At values of tan/3 > 52, the RGE evolution may 
break down due to a tendency towards a divergent bottom Yukawa coupling. 

The principal contributions to the total x 2 function for the (M A) tan (3) plane for scenario 
P2 are shown in Fig. 18. We see in panel (a) that the value of a M is very satisfactory in 
a band running across the plane from (M A ,tan(3) ~ (100 GeV, 15) to ~ (400 GeV, 50). In 
particular, large values of M A and small tan/3 are disfavoured. Panel (b) shows that the 
LEP lower limit on M h disfavours M A < 300 GeV also in this scenario. Small values of M A 
are also disfavoured by b — > sj, as seen in panel (c), and large values of M A and tan/3 are 
also disfavoured, but to a lesser extent. Panels (d) and (e) show that large values of tan/3 
and small values of M A are again disfavoured by both B s — > and B u — > rv T . Finally, 

panel (f) shows the combined EWPO and BPO x 2 function for scenario P2 throughout the 
(M A , tan /3) plane. We see that the best-fit point has M A ~ 340 GeV and tan/3 ~ 35. This 
is a good fit to both a M and b — > 57, as well as the M h and other BPO constraints. We note 
that this best-fit point has x 2 — 3.5, which is a noticeable improvement on the CMSSM 
fits discussed in the previous Section. We note that the A% 2 = 2.30 and 4.61 contours are 
somewhat more compact than in the case of scenario PI. 

For the parameter choice of P2 large values of M A are excluded because the right-handed 
stau becomes the LSP. This could be avoided by lowering the value of mi/ 2 outside the range 
in eq. (41), so as to recover a neutralino LSP. However, unless we drop my 2 substantially 
below our adopted range, the relic density will be too small due to LSP-stau coannihilations. 
Finally, we note that the hole around (M A , tan /3) ~ (600 GeV, 17) is due to the funnel. 
In this hole, the relic density is far too small to supply the preferred amount of cold dark 
matter. However, the hole could be filled if a larger range were chosen for my 2 - 

The principal contributions to the total x 2 function for the (M A: tan /3) plane for scenario 
P3 are shown in Fig. 19 10 . We see in panel (a) that the value of a M is satisfactory only for 
very large values of tan/3, almost independently of M A . In particular, values of tan/3 < 25 
are quite strongly disfavoured. Panel (b) shows that the LEP lower limit on M h disfavours 
very low values of tan (3 and M A < 150 GeV for tan (3 < 25. Small values of M A and tan (3 
are also disfavoured by b — > 57, as seen in panel (c), as also are large values of M A and 
tan (3. Panels (d) and (e) show that large values of tan (3 and small values of M A are again 

10 We do not display the x 2 values in the underdense slivers of the plane. 
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Figure 17: The most important contributions to the total x 2 value for the NUHM 
(Ma, tan (3) plane PI, due to (a) a M , (b) M h , (c) b — > 57, (d) B s — > and (e) B u — > ru T , 

and (f) the combined EWPO and BPO x 2 function. We use m t = 171.4 ±2.1 GeV and 
mb(mb) = 4.25 ±0.11 GeV, and mi/2 is adjusted continuously so as to yield the central value 
of the cold dark matter density indicated by WMAP and other observations for the central 
values of m t and ra&(rab). 

31 



55 
50 
45 
40 
35 

I 30 
25 

20 

15 

10 



3<i <5 
2 < X < 3 
1 < t < 2 

o < x z < 1 



••••••••••• 



ttt i nm tiss tt iisii issi a 



100 200 300 400 500 600 700 800 900 1000 



•••••••••••••A 

••••••••••••••• 



'1 1 1 1 1 1 

x'-K 



s<x ■ 
3<x'- 
2<x*. 



100 200 300 400 500 600 700 800 900 1000 



X : BR(b -> sy) 
li = 800, m n = 300 



••••••••••••••a******** 



••••••••••••••••••a 



••••••••••••• 



5 K i i i I i i i i I i i TfTTtTTttTTytTTttytTf T , I I 

100 200 300 400 500 600 700 800 900 1000 
M fl 



••••«••••• 



m 1 1 1 1 1 1 1 1 1 1 1 

X 2 : BR(B s ->nV) 
p = 800. m D = 300 
• 10<i s 



u 1 1 1 1 1 1 ttt ttttt ttmttmt tttttt mImmImmL- 

100 200 300 400 500 600 700 800 900 1000 



•••••• •••••••••••• 



••••••••••••• 

•••••••••••••« 

••••••••••a 



••••• ••••• 



1 1 1 1 1 1 1 1 ttt i mmm t Tttmttm t , , ■ i , , , , i ■ , , , i ■- 

100 200 300 400 500 600 700 800 900 1000 
M fl 



l t = 171.4, x m!n = 3.5 



• AX < 2 30 
X best fit 



100 200 300 400 500 600 700 800 900 1000 



Figure 18: The most important contributions to the total x 2 value for the NUHM 
(Ma, tan (3) plane P2, due to (a) a M , (b) M h , (c) b — > 57, (d) B s — > and (e) B u — > ru T , 

and (f) the combined EWPO and BPO x 2 fucntion. We use m t = 171.4 ±2.1 GeV and 
nib(mb) = 4.25 ±0.11 GeV, and mi/2 is adjusted continuously so as to yield the central value 
of the cold dark matter density indicated by WMAP and other observations for the central 
values of m t and m,b(mb). 
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disfavoured by both B s — > and B u — > rz/ r n . We note that AM Bs (not shown) also 

disfavours small Ma and large tan (3. Finally, panel (f) shows the combined EWPO and 
BPO x 2 function for scenario P3 throughout the (Ma, tan (3) plane. We see that the best-fit 
point has M A ~ 300 GeV and tan/5 ~ 35. This is not a very good fit to a M , but it is a 
good fit to b — > S7, M h and the other BPO constraints. We note that this best-fit point has 
X 2 = 7.4, which is not an improvement on the CMSSM fits discussed in the previous Section. 
The Ax 2 = 2.30 and 4.61 contours are somewhat looser than in the two previous scenarios, 
and extend to very large M A - 

The principal contributions to the total x 2 function for the (Ma, tan (3) plane for scenario 
P4 are shown in Fig. 20 12 . We see in panel (a) that a M favours a swathe with tan (3 ~ 15 to 
20, almost independently of M A - In particular, values of tan/3 > 25 and tan/3 < 5 are quite 
strongly disfavoured. Panel (b) shows that the LEP lower limit on M h disfavours values 
tan/3 < 15, the constraint becoming stronger for Ma < 200 GeV. A small band of values of 
Ma and tan/3 are favoured by b — > sj, as seen in panel (c), extending to tan/3 > 15 only for 
M A below the funnel at ~ 250 GeV. Panels (d) and (e) show the familiar feature that large 
values of tan/5 and small values of Ma are disfavoured by both B s — > fi + fi~ and B u — > tu t , 
and the same is true for AM Ba (not shown). Finally, panel (f) shows the combined EWPO 
and BPO x 2 function for scenario P4 throughout the (M A ,tan(3) plane. We see that the 
best-fit point has M A ~ 200 GeV and tan/3 ~ 20. This is not a very good fit to M h , but it 
is a good fit to a M , b — > 37, and the other BPO constraints. We note that this best-fit point 
has x 2 — 5-6, which is similar to the CMSSM fits discussed in the previous Section. We note 
also that the Ax 2 = 2.30 and 4.61 contours are particularly tight in this scenario, and rule 
out very large values of Ma and/or tan (3. 

There are some common features of these analyses for fixed (mi/2, too) and (mo, fx). For 
example, we find that relatively low values of M A ~ 200 to 400 GeV are consistently favoured. 
This is essentially because a M prefers moderately small values of TO1/2 which would, if left to 
themselves, create problems for b — > sj. However, this tension may be mitigated if Ma is 
correspondingly small, providing a cancellation in the supersymmetric contributions to the 
b — > 57 decay amplitude. We also note a consistent preference for relatively large values of 
tan (3 ~ 20 to 50, which is essentially due to the pressure exerted by the LEP lower limit on 
the Higgs mass. 

As discussed above, the LSP would constitute (most of) the cold dark matter across 

11 We note in panel (e) the appearance of a second, narrow favoured strip of parameter space. In this strip, 
the charged- Higgs contribution to the decay amplitude is not a small perturbation, but is <~ — 2x the W ± 
contribution! 

12 Again, we do not display the \ 2 values in the undcrdense slivers of the plane. 
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Figure 19: The most important contributions to the total x 2 value for the NUHM 
(Ma, tan /3) plane P3, due to (a) a M , (b) M hj (c) b — > 57, (d) 5 S — > and (e) B u — > rv T , 

and (f) the combined EWPO and BPO \ 2 function. We use m t = 171.4 ±2.1 GeV and 
mb(mb) = 4.25 ±0.11 GeV, and fx is adjusted continuously so as to yield the central value 
of the cold dark matter density indicated by WMAP and other observations for the central 
values of m t and 7715(7715). 
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Figure 20: The most important contributions to the total x 2 value for the NUHM 
(Ma, tan (3) plane P4, due to (a) a M , (b) M h , (c) b — > 57, (d) B s — > and (e) B u — > ru T , 

and (f) the combined EWPO and BPO \ 2 distribution. We use m t = 171 A ±2.1 GeV and 
mb(mb) = 4.25 ± 0.11 GeV, and fi is adjusted continuously so as to yield the central value 
of the cold dark matter density indicated by WMAP and other observations for the central 
values of mt and mb(mb)- 

35 



(most of) the NUHM parameter planes discussed above. Accordingly, for completeness we 
discuss the prospects for direct dark matter detection in different regions of the planes. In 
the cases of planes PI and P2, the direct scattering rate is generally below the CDMS upper 
limit [111], once one takes into account uncertainties in the strange-quark contributions to the 
spin-independent scattering matrix elements and in the local cold dark matter density. In the 
cases of plane P3, only in the region where Ma is small and tan f3 is high does the dark matter 
scattering rate approach the CDMS upper limit. However, there is a potential conflict with 
the preliminary XENON10 results [112] if the strange-quark contribution and/or the local 
relic density is large. A similar situation arises in plane P4 for small values of Ma, almost 
independent of tan j3. In absence of sufficient understanding of the systematic uncertainties 
in the strange-quark contribution and the local cold dark matter density, we do not attempt 
to include the direct dark matter searches in the overall x 2 function 13 . 

The survey of NUHM parameter space made in this Section has not been exhaustive, in 
particular we have restricted our attention to planes with A = 0. Nevertheless, the values of 
X 2 found at the best-fit points in the various (M^,tan/3) planes are quite acceptable: planes 
PI, P2, P3 and P4 have x 2 — 7.1, 3.5, 7.4 and 5.6, respectively, in fits to 9 observables with 
2 free parameters in each case. It should be stressed, however, that only the P2 plane has 
a minimum of x 2 noticeably lower than that for the CMSSM fits with tan (3 = 10, which 
occurs when /x = 800 GeV, m = 300 GeV, M A ~ 340 GeV and tan (3 = 36 and, moreover, 
at the point with the minimum value of x 2 , the relic neutralino density is somewhat higher 
than the WMAP-compatible range. One might expect a greater reduction in x 2 in a full 
study of the NUHM, in view of its two additional parameters compared with the CMSSM. 
Accordingly, we have made a preliminary study whether the quality of the NUHM fit could 
be improved significantly by varying A , assuming the same values of mo, Ma and tan/3 as 
at the best-fit point in the P2 plane, but choosing different values of /x and m 1 /2- We have 
investigated the possibilities (n,m 1/2 ) = (800,368) [LI], (800,448) [L2] and (680,448) GeV 
[L3], respectively, and varied A between ±1000 GeV. Fig. 21 shows the values of x 2 along 
the lines LI, L2 and L3. We see that the greatest improvement in x 2 compared to the 
(Ma, tan P) planes shown previously are by ~ 0.3 only. Interestingly, the minimal values of 
X 2 are found for small values of A ~ 0. Undoubtedly some further reduction in x 2 could 
be found in a more complete study, but it seems that the extra degrees of freedom in the 
NUHM are not crucial for the overall quality of the fit. 

13 For completeness, we note that along the WMAP strips in the CMSSM the direct dark matter scattering 
rate is always comfortably below the CDMS upper limit. 



36 



14 



12 



~i — i i | i i i | i i i | — i i i | i i i | i i i — prn — i | i i i | i i — r~p~i — r~r 

NUHM, M A = 340 GeV,tanp = 36, m t = 171.4 GeV - 

m 1/2 = 368 GeV, \i = 800 GeV 

- m 1/2 = 448 GeV, ^ = 800 GeV, CDM too large 

m 1/2 = 448 GeV, ^ = 680 GeV 




J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L_l_ 

000 -800 -600 -400 -200 200 400 600 800 1000 

A [GeV] 



Figure 21: The dependence on Aq of the % 2 function along lines with mo = 300 GeV, Ma = 
340 GeV, tan/? = 36 and (n,m 1/2 ) = (800,368) [Ll], (800,448) [L2] and (680,448) GeV 
[L3], respectively. 

5 Conclusions 

We have analyzed previously the regions of the CMSSM parameter space preferred by the 
EWPO [8,9], and found a tendency to prefer regions on the WMAP coannihilation strips with 
relatively low values of m\/2- These points were favoured, in particular, by the measurements 
of a M and My/- Both these tendencies have now been reinforced, with the interpretation of 

based on the use of e + e~ data to estimate the SM contribution gaining ground, and the 
small decrease in m t and the slight increase in Mw tending to favour a contribution to the 
latter EWPO from some physics beyond the SM. 

Previously, we incorporated just a single BPO into our global analysis [8,9], namely 
b — > sj. Recently, data on B s — > B u — > tv t and AMb s have also become available 

and now impinge significantly on the CMSSM parameter space. In this paper, for the first 
time, we have incorporated all these BPO into a global analysis. 

We have found a good x 2 /d-o.f. for this global fit. However, it is clear that there is 
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a slight tension between the relatively low values of mi/ 2 favoured by the EWPO and the 
absence of any corroborating indication from the BPO. Nevertheless, the global x 2 analysis 
favours the appearance of relatively light sparticles that should be 'easy' to see at the LHC 
and may offer good prospects also for the ILC. 

As we also have discussed here explicitly, for the first time, the global analysis strongly 
favours values of M h only slightly above the lower limit established by LEP. Indeed, we find 
that values of M h < 120 GeV are preferred, while values above 123 GeV cannot be reached 
in the CMSSM (for the current m t value). 

Another new step in this paper, motivated by the slight tension between the EWPO and 
the BPO, has been to explore the parameter space of the NUHM. We have displayed for 
the first time (Ma, tan (3) planes within the NUHM over which the WMAP constraint on 
the cold dark matter density is generically respected. We have then shown the interplay of 
the various EWPO and BPO in such planes. We find that, for fixed (m 1 / 2 ,m ) or (m ,/i), 
relatively low values of Ma ~ 200 to 400 GeV are favoured, as are relatively large values of 
tan (3 ~ 20 to 50. It is possible to find in this way NUHM points that have lower x 2 than 
those possible in the CMSSM. 

In the future, it will be necessary to follow closely the evolutions of both the EWPO 
and the BPO: improvements in the measurements of both M w and m t are expected, and 
the interpretation of a M may become clearer. We also expect significant improvements in 
the measurements of B s — > and B u — > tu t , and possibly in the interpretations of 

b — > S7 and AM Ba - The question will be whether the present slight tension between the 
EWPO and the BPO within the CMSSM will strengthen or relax, and a more detailed 
and systematic exploration of the NUHM parameter space will certainly be desirable. The 
implications of the EWPO and the BPO for the supersymmetric parameter space will surely 
be an interesting and continuing saga. 
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